Detection of high-field superconducting phase in CeCoIn 5 with magnetic susceptibility 
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We measure the ac susceptibility of single-crystal CeCoIns in dc field parallel to the c axis and 
find further evidence for a high-field phase transition within the superconducting phase in this 
orientation. We apply up to 2.3 kbar uniaxial pressure along the c axis and discuss the pressure 
dependence of the high-field phase. We also report the behavior of H C 2 under uniaxial pressure for 
field along c. 



The heavy-fermion superconductor CeCoIn5 and its 
isostructural cousins Ce(Ir,Rh)In 5 and Pu(Co,Rh)Ga 5 
exhibit several intriguing features. Collectively these 
compounds have ground states that are superconducting, 
antifcrromagnetic, or both. The materials have signifi- 
cant two-dimensional character, suggested by their lay- 
ered crystal structure and confirmed through electronic 
structure calculations. Furthermore, within each family 
T c varies linearly with the ratio of lattice constants c/a, 
illustrating how closely geometry and superconductivity 
can be linked [l|. With a c-axis magnetic field CeCoIn5 
also has a quantum critical point near the upper criti- 
cal field, which can be shifted into the superconducting 
phase with hydrostatic pressure @- 

Yet another unusual feature appears in the high-field, 
low-temperature portion of the superconducting regime. 
This has been studied mainly for magnetic fields applied 
in the basal plane, where H C 2 is about 11.8 Tesla and a 
signature within the superconducting phase is seen near 
10 Tesla with magnetization Q, specific heat Q, neu- 
tron diffraction and NMR Q measurements. The 
additional signature moves to higher fields as tempera- 
ture increases, intersecting the H C 2(T) curve near 300 
mK. These data are of particular interest since CeCoIns 
is a good candidate to support the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) superconducting phase described 
almost 50 years ago @, H|- This state is marked by 
a non-zero quasiparticle center-of-mass momentum and 
a spatially periodic superconducting order parameter. 
CeCoIn.5 is in the extremely clean limit and has strong 
Pauli limiting of its critical field. In addition, the high- 
field phase (HFP) is accompanied by a change of the 
normal-superconducting transition from second-order at 
low fields to first-order at high fields, in agreement with 
FFLO theory 0, i, \M- While the exact nature of the 
phase and its identification with FFLO remain in ques- 
tion 0, 0] , the experimental evidence for the HFP with 
in-plane field is clear. 

Much less effort has been devoted to CeCoIns with 
magnetic field along the c axis. A few measurements 
have indicated a high-field phase for this orientation 
too, as seen both by signatures within the supercon- 
ducting phase [LlTtl3l | and by a change of the normal- 



superconducting transition from second-order to first- 
order similar to that observed for in-plane field. How- 
ever, other experiments have found no signatures 0,0], 
leaving even the experimental situation unclear. In this 
report we provide further evidence for the HFP phase in 
CeCoIns for H\\c, along with the uniaxial pressure de- 
pendence of the HFP boundary up to 2.32 kbar. 

We measured three samples oriented with uniaxial 
pressure, dc field, and ac field parallel to the c axis. The 
pressure is applied with a bellows setup activated with 
helium gas from room temperature. The samples had 
mass of 3.67 mg, 0.67 mg and 0.92 mg with area 2.78 
xlO -6 m 2 , 5.16xl0~ 7 m 2 , and 6.97 xl0~ 7 m 2 respec- 
tively. The experimental details, including sample prepa- 
ration and orientation, are identical to those we reported 
previously with the exception of the additional dc 
field used in this study. 
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FIG. 1: Onset of superconductivity for different dc fields and 
at fixed pressure of 1.82 kbar, with overlaid fits. Legend in- 
dicates the applied field in Tesla for each curve. 



2 



For four runs in the pressure cell at nominally zero 
pressure and zero field, we found an average T c = 2.29 
K. In each case the transition temperature measured in 
the pressure cell was within 24 mK of the value outside 
the pressure cell. In some cases the temperature shift 
can indicate the initial pressure applied during the ex- 
perimental setup and cooling process, typically about 0.3 
kbar [14l. |15|. However, T c in CeCoIns is not very sensi- 
tive to c-axis pressure, shifting by less than 17 mK/kbar 
[lH liH ] , so it gives little information on the pressure off- 
set. Hence the values for applied pressure used in this 
paper do not include any offset for an initial pressure. 

Figure Q] presents susceptibility curves near the onset 
of superconductivity for several applied fields, all at a 
pressure of 1.82 kbar. As indicated by the black lines in 
the figure, we determine T c by assuming the susceptibil- 
ity is constant in the normal state and changes linearly 
just below T c . As usual for superconductors, the onset T c 
shifts to lower temperature as the applied field increases. 
In addition, the transition decreases in magnitude and 
broadens in temperature. The magnitude change is due 
to flux lines threading the sample, which lead to incom- 
plete diamagnetism. The broadening is governed by the 
temperature-dependence of the vortex pinning. As tem- 
perature falls and the vortices are more strongly pinned, 
the behavior becomes increasingly diamagnetic. We now 
turn to temperature sweeps below 400 mK and at fields 
above 4.5 Tesla. In this regime the susceptibility has not 
saturated but retains a temperature-dependence which 
decreases with increasing field. 

We confirmed that changing the pressure causes no sys- 
tematic shift in either \s or \ n . Nor does field shift Xn- 
Together these indicate that only \s is affected and only 
by tuning applied field. This corresponds to a decrease 
in superconducting transition magnitude with increasing 
field. 

Susceptibility at much lower temperatures appears in 
Figure [5] The data shown were measured at 1.42 kbar 
uniaxial pressure. When we increase the dc field the su- 
perconducting susceptibility Xs decreases in magnitude. 
We track this shift as a function of field by picking out \ 
at a fixed temperature from each of these dat sets. Fig- 
ure [3] shows some of the resulting x vs H curves. We 
repeat this procedure for each of the four temperatures 
indicated by the dashed lines of Figure [5J In each case 
we smooth the data by averaging % in a 40 mK range 
centered at the desired temperature. 

The upper frame of Figure |3] shows data at 300 mK for 
several different pressures. Each data set appears to have 
linear regimes, but also a kink where the slope changes 
abruptly. At each pressure we fit a function of exactly 
this form: two linear portions, with a kink at (Hk,Xk) 
where the slope changes. We use four free parameters: 
Hk, Xk, and the slope on each side. We vary Hj~ and 
Xk manually, in steps of 0.01 Tesla and 0.00075, respec- 
tively. For each kink location, we solve for the slopes 
which best fit the data above and below ■ We then se- 
lect the kink location which minimizes the least squares 
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FIG. 2: Susceptibility deep in the superconducting phase at 
1.42 kbar. The applied field for each curve is indicated to 
the right. The vertical dashed lines indicate the temperatures 
that we use to set up the field vs temperature curves as shown 
in Figure [3] 



error. The lines overlaid in the center plot are the result 
of this minimization. The lower plot indicates the loca- 
tion of Hk and also the critical field H C 2 at 300 mK for 
each pressure. The two decrease in concert with increas- 
ing pressure, showing that most of the kink's shift with 
pressure comes from the reduction of H c2 . 

The main portion of Figure @] shows H C 2 (T) at ambi- 
ent pressure, with also plotted. As seen in the lower 
part of Figure [31 c-axis pressure reduces H C 2 at low tem- 
peratures. At 500 mK, the change in H C 2(T — 0) is 0.04 
T/kbar. Hydrostatic pressure has a somewhat larger ef- 
fect, between 0.06 and 0.08 T/kbar j|,[I3|- Interestingly, 
although increasing hydrostatic pressure also pushes the 
high-temperature, low-field portion of the H C 2(T) dome 
up to higher T c , this shift is almost absent with uniaxial 
pressure. We recently reported that the zero-field T c has 
little dependence on c-axis pressure, rising only about 20 
mK to a maximum near 2 kbar, so it is not surprising 
that c-axis pressure also has little effect on the phase 
boundary at very low fields. 

The inset of Figure 2] expands the high-field, low- 
temperature region. The solid curve is H C 2(P = 0), nor- 
malized to 1 at T = 0. The normalized H C 2 curves at 
other pressures are nearly identical and are omitted from 
the graph. The symbols show the normalized kink field, 
H k (T, P)/H c2 (0, P). The kink agrees well with the pro- 
posed HFP boundary, which appears between 0.955 and 
0.96 in other measurements [a EI EBIlil- Several tran- 
sitions of the vortex lattice have also been observed [l8[ , 
but the highest of these is near 0.85, much lower than the 
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FIG. 3: Top: Susceptibility at 300 mK with fits overlaid. See 
text for details. All data except zero pressure are from the 
same sample and are shifted vertically for clarity. Bottom: 
Plot of the location of the kinks Hk obtained from the fits 
(filled circles) , with H C 2 at 300 mK also shown (open squares) . 
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FIG. 4: Main graph presents measurements of H C 2(T) at zero 
pressure (solid circles) and a fit to these points. Shaded tri- 
angles indicated the kink locations Hk(T). Inset shows an 
expanded view in the region of Hk near the top of the dome. 
The curve H C 2(T) is normalized to 1 at T — 0, and Hk is 
shown for several applied pressures. 



kink. Under uniaxial pressure both H c2 an d H k decrease, 
by roughly the same factor. Hk is also approximately in- 
dependent of temperature and forms a linear boundary 
well below H c2 . 

To complement the temperature sweeps described 
above we also performed field sweeps at fixed tempera- 



ture and pressure. One might expect a kink as in Figure 
[3] to appear directly in this measurement. However, in 
fact we find no evidence of a kink. At all fields below 
H C 2, the slope of x vs H in the field sweeps is compara- 
ble to its high-field value in the temperature-sweep data. 
The low-field slope in the temperature sweeps is notably 
larger. The sensitivity to the measurement history could 
explain why some of the previous experiments have found 
no sign of a transition to the HFP for c-axis field. 

One possible source of the kink is influence from the 
AFM QCP, which at ambient pressure occurs at a field 
very close to H c2 0. However, under hydrostatic pres- 
sure the QCP field decreases about five times as fast as 
H C 2, so that the QCP moves deep inside the supercon- 
ducting phase. While its behavior has not been tracked 
under uniaxial pressure, there is no reason to believe that 
H C 2 and Hqcp would remain close together. Since we 
find that Hk changes with uniaxial pressure only about 
1.5 times as fast as H C 2, the kink is probably not an 
indication of the QCP. 

Another interpretation is in terms of thermally acti- 
vated flux flow [20]. Our temperature sweeps are nearly 
field-cooled measurements. Although in practice we of- 
ten do not exceed T c , we find no apparent difference in 
the low-temperature susceptibility between these temper- 
ature sweeps and those when the sample does begin in 
the normal state. The larger low-temperature slope for 
temperature sweeps indicates that the vortex pinning is 
stronger after field-cooling. A plausible explanation is 
that at higher temperatures the vortices fluctuate enough 
to find the most favorable pin sites, which they cannot 
do when the field changes at low temperature. In the 
HFP, there is no such annealing effect, and the pinning 
strength always has the smaller value. 

A change in pinning strength could arise from antifer- 
romagnetic order in the region above Hk, which could 
interact with the superconductivity. More generally, any 
change in the order parameter can reasonably be ex- 
pected to affect the vortex pinning. Several reports have 
been made along with theoretical descriptions suggest- 
ing this phase is the FFLO phase (see reference [19| for 
an overview). In the FFLO state periodic planar nodes 
appear perpendicular to the flux lines, leading to a seg- 
mentation of the vortices into pieces of length A = 2n/q, 
where q is the quasiparticle wave vector [l9| which fac- 
tors into the order parameter A oc sinq • r. The pieces 
are relatively free of each other and hence better able 
than conventional vortices to position themselves at pin 
centers. Hence pinning vorces on the flux lines should in- 
crease in the FFLO phase. If the nodal planes reflect the 
crystal structure, then uniaxial pressure parallel to the 
flux lines would compress the segmentation length A and 
change the periodicity of the order parameter. Already 
strongly Pauli-limited, CeCoIns under uniaxial pressure 
becomes even more so than with hydrostatic pressure, 
possibly due to increased hybridization between the Ce- 
In layers. The way that Hk tracks H c2 suggests a possible 
connection with the electron spin population. Theoret- 
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ical models suggest redistribution of spin states in the 
vortex cores upon entering the FFLO state [2l[ . Namely, 
at the intersection of the nodal plane and the vortex core, 
excess spin states are emptied. 



Ultrasound measurements of the HFP boundary [22 1 
for in-plane field can also be interpreted in terms of vor- 
tex pinning. The vortices affect the ultrasound velocity 
more strongly in the HFP, which is consistent with in- 
creased pinning; one explanation is segmentation of the 
vortices into short, somewhat independent pieces that 
can better take advantage of low-density pin sites. Such 
segmentation could arise from the planar order parameter 
nodes of an FFLO state. By contrast, our work suggests 
lower pinning in the HFP. This may be evidence that 
the HFP phases for the two field orientations are not 
the same. However, we note that the two measurements 
are not necessarily contradictory. The ultrasound signal 
comes from the most strongly pinned vortices, while our 
measurements probe the vortices that are most free to 
move. In principle the HFP could support a wider range 
of pinning strengths that produces both effects. 

In conclusion, we have measured ac susceptibility re- 



sponse of CeCoIns in an applied field, as well as under 
uniaxial pressure. We find evidence for a phase boundary 
in agreement with previous indications of the transition 
to the HFP at zero pressure. We report a shift in this 
boundary with uniaxial pressure which roughly tracks 
H C 2, decreasing in field at a rate about 50% faster than 
the decrease of H c2 . For H c2 itself, pressure depresses 
the critical field at low temperatures but has little effect 
near T c . This differs from hydrostatic measurements but 
is consistent with the minimal effect of uniaxial pressure 
on T c itself. The location and pressure dependence of the 
fife boundary provide further evidence for the HFP phase 
for H\\c in CeCoIns. An unanswered question is why the 
iJfc boundary shows up with temperature sweeps and not 
field sweeps. 
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